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Abstract 

In this paper, we have investigated the entanglement between two dipole coupled 
two-level atoms. The model, in which only one atom is trapped in an lossless 
cavity and interacts with single-mode thermal field, and the other one can be 
spatially moved freely outside the cavity has been carried out. We have consid¬ 
ered the effect of the atomic coherence on the entanglement behavior. We have 
shown that a thermal field might cause high entanglement between the atoms 
both for coherent and incoherent initial atomic states only for small values of 
the cavity mean photon number. In the considered model the atoms would get 
entangled even when both atoms are initially in the excited state. We have also 
derived that the degree of entanglement is weakly dependent on the strength of 
dipole-dipole interaction for coherent initial states. 

Keywords: Entanglement, Cavity QED, Negativity, Thermal field, Atomic 
coherence 


1. Introduction 

Entanglement between separate quantum systems is one of the key problem 
in quantum mechanics. It plays a central role in quantum information, quan¬ 
tum computation and communication, and quantum cryptography [l|. Sev¬ 
eral methods of creating entanglement have been proposed involving trapped 
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and cooled ions or neutral atoms in cavities, superconducting circuits, spins 
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in solids etc. [2f. In order to function optimally various applications require 
maximally entangled states. Because of decoherence, which is generally related 
to noise, there is great difficulty in generating and keeping the integrity of a 
pure entangled states. Although the interaction between the environment and 
quantum systems can lead to decoherence, it may also be associated with the 
formation of non-classical effects such as entanglement [sj. Thus, understand¬ 
ing and investigating entanglement of mixed states becomes one of the actual 
problem of quantum information. Recently, Bose et al. [3 have shown that 
entanglement can always arise in the interaction of an arbitrary large system in 
any mixed state with a single qubit in a pure state, and illustrated this using 
the Jaynes-Cummings interaction of a two-level atom in a pure state with a 
field in a thermal state at an arbitrary high temperature. Kim et al. 0 have 
investigated the atom-atom entanglement in the system of two identical two- 
level atoms with one-photon transition induced by a single-mode thermal field. 
They showed that a chaotic field with minimal information can entangled atoms 
which were prepared initially in a separable state. Zhang directly generalizes 
Kim’s study to the case when the atoms are slightly detuned from the thermal 
field, and study how the detuning would affect atom-atom entanglement 0. 
Zhou et al. [§] have considered the same problem for nonidentical atoms with 
different couplings. The entanglement between two identical two-level atoms 
through nonlinear two-photon interaction with one-mode thermal field has been 


studied by Zhou et al. [7]. They showed that atom-atom entanglement in¬ 
duced by nonlinear interaction is larger than that induced by linear interaction. 
In fl has discovered that two atoms can be entangled also through nonlinear 
nondegenerate two-photon interaction with two-mode thermal field. The influ¬ 
ence of dipole-dipole interaction on entanglement between two cubits induced 

iqi-H- 


by one-mode and two-mode thermal field has been investigated in 


The problem of creating or controlling the atomic entanglement is greatly 
related to the atomic coherence of population between different levels. The 
authors of papers 14j-[18] have shown that the entanglement between two atoms 
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induced by one-mode or two-mode thermal field can be manipulated by changing 
the initial parameters of the atoms, such as the superposition coefficients and 
the relative phases of the initial atomic coherent state and the mean photon 
number of the cavity field. They have also discovered that entanglement may be 
greatly enhanced due to dipole-dipole interaction in the presence of the atomic 
coherence 171. 

While, . has hee„ m e„«i„„ed ip Q ,he prac.ica, applied ip ,„a„- 
turn information processing require engineering entangled atoms so this expects 
operable atoms which can be moved to distance without losing of information. 
Recently many schemes have been proposed to realize the engineering entangled 


atoms 


20j-[25]. Guo and coauthors have proposed a simple scheme to realize 


an easily engineered two-atom entangled state. The advantage of this scheme 
is only one atom is trapped in a cavity, and the other one can be spatially 
moved freely outside the cavity. But the authors have investigated the atom- 
atom entanglement only for vacuum initial cavity field. In this paper we study 
the entanglement dynamics for model supposed in [19{] for thermal cavity field 
taking into account the initial atomic coherence. 


2. The model 

We consider two identical two-level atoms and one-mode quantum electro¬ 
magnetic cavity field. The first atom is trapped in a lossless microcavity and 
resonantly interacts with the cavity field of the frequency w. The second atom 
lies beside the first atom out of the cavity. We assume that the distance between 
atoms can compare with a wavelength on working transition. In this case the 
dipole-dipole interaction should be included. The Hamiltonian of this system 
can be written as 

H = (l/2)hui(al + U 2 ) + hiua + a + hg(af a + a+cr-f) + hJ(af ), (1) 

where (l/2)of is the inversion operator for the zth atom (z = 1,2), af = 
|) and of - = |— }ii(+| are the transition operators between the excited 
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|+)i and the ground |—)i states in the zth atom, a + and a are the creation 
and the annihilation operators of photons of the cavity mode, g is the coupling 
constant between atom and the cavity field and J is the coupling constant of 
the dipole interaction between the atoms and |+) and |—) are the excited and 
the ground states of a single two-level atom. The two-atom wave function can 
be expressed as a combination of state vectors of the form | vi,v%) = |iq||i > 2 ), 
where Vi , ug = +, —. 

We consider that the initial states of atoms are the coherent superposition 
of the two levels, that is, 

|Ti(0)) = cos0i|+)i + e* Vl sin6>i|—}i, 

|’MO)) = cos0 2 |+)i + e* V2 sin1—) 2 - (2) 

Here 9\ and 9 2 denote the amplitudes of the polarized atoms, and and </j 2 
are relative phases of two atoms, respectively. So the initial density matrix for 
two-atom system can be written as 

Pa( 0) = |M0)}|M0))(M0)|(M0)|. 

The initial cavity mode state are assumed to be the thermal one-mode state 

Pf( 0) = '^2pu\n)(n\. 

n 

The weight functions are 

n n 

Pn ~ (1 + n)” +1 ’ 

where h, is the mean photon number in the zth cavity mode, 

fh = (exp [hiOi/k B T} - l] -1 , 

ks is the Boltzmann constant and T is the equilibrium cavity temperature. 
The initial density matrix of the whole system is 

P( 0) = Pa(0) Pf (0) = ^p„|^ 1 (0))|v|/ 2 (0))(vl> 1 (0)|(^ 2 (0)||n)(n|, (3) 


4 



Before considering the dynamics of the system for thermal initial cavity 
field, it is straightforward to first study the case when the trapped two-level 
atom interacts with Fock state. Suppose that the excitation number of the 
atom-field system is n ( n > 0) the evolution of the system is confined in the 
subspace | —, —, n + 2), |+, —, n + 1), | — n + 1), |+, +, ra). On this basis, the 
eigenfunctions of the Hamiltonian (1) can be written as [191 ] 

l^ira) = ^in (-^Qlro ] > —, tl + 2) + Xi 2n \ + , —, 71 + 1) + 


+^i 3 n|~) +, n + 1) + Xj 4 n |+, +, n }) (* — 1,2,3,4), 


where 

£in = 1/\/|^il rl | 2 + |Xj2n| 2 + |A"i3„| 2 + |A"i4n| 2 


and 

x. iln = \, x i2n = {-i) i+1 -^=, 

V 71 2 

_ Si - (re + 2) 1 ^ +1 SnjSn — (re + 2) — a 2 ) 

(\/n + 2)a ’ ?4 " (\/n + l\/?z + 2)a 

Here ;S' n = for odd i and S n = B n for even i. The corresponding eigenvalue 


are 


E\ n /h — (re + l)u> + A n , E 2n /h — (re + l)w — H n , 

E 3n /h = (re + 1 )u> + B n , E 4n /h = (re + l)w - B n . 


Here 

An = V Wn + V n /2, B n = y/W n ~ V n /2 

and 

W n =4n + 6 + 2a 2 , F n = 2\/4(n + l)a 2 + (a 2 + 1), 
where a = J/g. 

To derive the full dynamics of our model one can consider also the basis states 
| —, —, 1), |+, —, 0), | —, + , 0). In this basis the eigenfunctions and eigenvalues of 
the Hamiltonian (1) are 


ki) = (a 2 /9)[| —, -, 1} - (l/a)| —, +, 0)], Ei = 0; 


5 













\<P2) = (i/v^)[(i/n)|-, i) + 1 +, -,o> + (a/n)|-, +,o)], e 2 /h = fi ; 

1^3) = (1/a/ 2 ) [—(i/fi)|—, 1) + 1 +, 0 ) - (a/n)|-, +, 0 )], E 2 /n = -n, 

where fl = \/l + a 2 . 

At last, the Hamiltonian (1) has one more eigenfunction 

<Po = | —, 0> 

which corresponds to energy Eq = —hui. 

Assume that the whole system is initially in the state |+, +,n) (n > 0), then 
at time t, the whole system will evolve to 

|'I'(t)} = Zu,n\ — , ~j tl+2)+^21,n| + > ^+l) + ^31,n| — t +) tl+l)+^41,n|+) + 5 n). 

(4) 

Here 

Zu,n = e~ lElnt ^ h Cm Y 41n x lln + e-' E *'V R & n Y 42 n X 21n + 
+e -iE 3 n t/H Y 43n x 31n + e ~ lEint ^ h Un Y 44n X 41n , 

Z 21 ,n = e~ lElnt / n Cm Y 41n X 12n + & n Y 42n X 22n + 

+e -iE 3 n t/h £ 3n Y 43n Xs2n + e ^E in t/h ^ Y 44n x 42n , 

Z 31 r n = e lElnt / h £ ln Y 41n A'i3n + e lE2nt / h £ 2n Y 42n X 23 n + 
+e -,B3nt/R Y 43n X 33n + e ~ lEint / n i 4n Y 44n X 43n , 

Z 4 l, n = e- lE ^ n Cm Y 4ln X 14n + e~ tE ^ h & n Y 42n X 24n + 
+ e -'E^t/n £ 3n Y 43n x 34n + e~ lEint / n £ 4 „ Y 44n X 44n , 
where Y ijn = { fn X* in . 

Similarly, when the whole system is initially in the state |+, —, n+1) (n > 0), 
then at time t, the whole system will evolve to 

|4'(i)} = Z\ 2 n | , —, n+2) + Z 22jn \ +, — > ?l+l) + -Z’32,n| — > + j tl + 1) + Z 42jH \ + j +j Tl). 

(5) 

Here 

Zi2, n = e lElnt / h Ci n Y 2 i n X 44n + e lE2nt / h £ 2 „ Y 22 n X 2 i n + 
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+e~' E ^ n b„. Y 23n X 31n + e - lE ^ h bn Y 24n X 41n , 

Z22, n = e lElnt / h £ ln Y 21n X\ 2n + e lE2nt / h \ £ 2n y 22n x 22n + 

+e -iE 3n t/h £ 3n y 23n x 32n + e -rE in tm ^ y 24n x 42n> 

Z 32 , n = e ~ lElnt/ ' H bn Y 21n X 13n + e - lE ^' n bn F 22 „ X 23n + 

+ e ~ lEsnt,h bn Y 23n X 33n + e -^t/n i4n Y 24n x 43n , 

Z 4 2,n = e~ lElnt / h bn Y 21n X 44n + e~ lE2nt ^ h bn Y 22n X 24n + 

+e~ lE3 " t/n bn Y 23n X 34n + e~ lE ^ t/n Un Y 24n X 44n . 

If the initial state of considered system is |+, —, 0), the time dependent wave 
function takes the form 

I *(*)> = Zn\-,-, 1) + Z 22 |+, 0) + Z 32 |- +, 0), (6) 

where 

% %Oi 

Z 42 = — — sin fit, Z 22 = cos fit, Z 32 = — — sin fit. 

For initial state | —, +, n + 1) (n > 0) the time-dependent wave function is 

|'F(f)} = Zi 3tn \ —, —, n+ 2 )+Z 23 ^ n \+, —, ?i+ 1)+Z33 jTI | —, +, n+l)+Z43 jri |+, +, n). 

(7) 

Here 

7 _ 0 —lEint/h t v V I „ — lE2nt/h f y y I 

^13,n — e Un'31n A HnTe <52n 2 32n^21n + 

+e~ lE3 " t/h bn Y 33n X 31n + e~ lE ^ h bn Y 34n X 41n , 

Z 23 ,n = e~ lElnt ^ h bn Y 34n X 42n + e _ ' E2n */ r, l bn Y 32n X 22n -\- 
+e ^E 3n t/n ^^ y 33h x 32n + e -' E ^l n bn Y 34n X 42n , 

Z 3 2,n = e~ >Elnt ^ H bn Y 31n x 13n + e-'&»*/* bn Y 32n X 23n + 

_l_ e lE 3 ^ t / n b n Y 33n X 33n +e lEint / n bn Y 34n X 43n , 
z 42 ,n = e~ lEln *^ h bn Y 34n X 44n + e ~ lE2nt / h bn Y 32n X 24n + 
+ e - tE ^l H bn Y 33n x 34n + e~ lE ^/ h bn Y 34n X 44n . 
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If the initial state of considered system is | —, +, 0), the time dependent wave 
function of the whole system takes the form 


l'k(f)} — Z i 3 | — , —, 1) + Z 23 \+, — j 0) + Z 33 \ — , +, 0), (8) 


where 


Z \ 3 = (cos Sit - 1), Z 23 = - ^ sin fit, Z 33 = -^-(1 +a 2 cos fit). 


For initial state | —, —, n + 2) (n > 0) we have at time t that 


|’F(f)} — Z i4 >n | —, —, n+2)+Z24 jn |+, —, ?r+l)+^ 34 >n | —, +, n-\-l)-\-Z 44}n \+, +, n). 

(9) 

Here 

Zl4,n = e~ lElnt / h Y 41n X lln + e- lE2 ^ n Y 42n X 21n + 

+e ^E 3 n t/n ^ y 43n Xain + e -zE in t/H Un y 44n X4inj 

Z 24 , n = e lElnt / h £ ln Y 41n X 42n + e £ 2n Y 42n X 22n + 

+e ^E 3 n t/n £ 3n y 43n X32n + e - tE 4 n t/h Un y 44n X42nj 

z 4 2 ,n = e~ lEln *^ h £1 n Y 41n X 43n + e ~ lE2nt / h £ 2n Y 42n X 23n + 

+ e ~ lE3nt/h ?3n Y 43n X 33n + e- lE ^/ R Un Y 44n X 43n , 

Z 44 ,n = e~ tElnt / n Y 41n X 14n + e~ tE2 ^ n & n Y 42n X 24n + 

+e -iE 3 n t/n ^ 3n y 43n Xa4n + e —iE 4 n t/fr Un y 44n X44n _ 

If the initial state of considered system is | —, —, 1), the time dependent wave 
function of the whole system takes the form 

|'F(f)} = Zi 4 \ — ,—, 1) + Z 24 \+, —, 0) + Z 34 \ — , +, 0), (10) 


where 

Z i 3 = {a 2 + cos fit), Z 23 = sin fit, Z 33 = --^(1 - cosClt). 

SZ S2 12- 

At last, if the initial state is | —,0), the time dependent wave function will 

evolve to 


|'F(f))=e-^ t |- - 0). 


( 11 ) 



Now we go back to the theme of this paper. If the initial state for two atoms 
is (2), then using the equations (4)-(ll) one can obtain the density operator for 
the whole system. Taking a partial trace over the heat bath variable one can 
obtain that the reduced atomic density operator in the two-atom basis evolves 
to 


PA{t) = 



' P li 

P12 

P 13 

P14 


Pl2 

P22 

P23 

P24 


Pis 

P23 

Pss 

P34 


y Pl4 

P24 

P34 

p44 

are given 

in Appendix. 


( 11 ) 


For two-qubit system described by the density operator pA{t ), a measure of 
entanglement or negativity can be defined in terms of the negative ei gen values 


/q of partial transpose of the reduced atomic density matrix p 1 ^ 26}, 27]. The 


partial transpose of the reduced atomic density matrix (11) is 

\ 


Pa M = 


The negativity is 


^ Pll 

Pl2 

Pl3 

P23 

Pl2 

P22 

P*14 

P24 

Pl3 

Pl4 

P33 

P34 

\ P23 

P24 

P34 

P44 

£ = - 

■2E 

Pi ■ 



( 12 ) 


When e = 0 two qubits are separable and e > 0 means the atom-atom entan¬ 
glement. The case e = 1 indicates maximum entanglement. 

The results of calculations of entanglement parameter (12) for initial atomic 
states (2) are shown in Fig. 1-3. For all figures we put (pi = = 0. 
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3. Results 


In Fig. 1 we plot the negativity as a function of gt for coherent and three 
incoherent atomic states described by formula (2). The curves corresponds 
to a fixed value of dipole strength a = 0.5 and weak thermal two-mode field 
with fi = 0.01. From these figure we can find that, firstly, entanglement can 
be induced by thermal field for all initial atomic state. The more significant 
result is that the atomic entanglement may be induced by thermal field when 
both the atoms are prepared in their excited states. By contrast, the atomic 
entanglement does not induced by thermal field when both atoms are trapped in 
cavity and prepared in their excited states Q . Secondly, the maximum degree of 
entanglement is slightly enhanced owing to the atomic coherence comparing the 
curves in Fig. 1(c) and Fig. 1(d). With increasing of the mean photon number 
(see Figs.2(a) and Fg2(b)) the value of atom-atom negativity decreases. But for 
coherent states this decreasing is much sharper. Such entanglement behavior 
essentially differs from that for model when both atoms are trapped in cavity 
and interact with a thermal field. In the last case jji| the entanglement is 
greatly enhanced due to the initial atomic coherence for thermal cavity with 
large mean photon numbers. Fig. 3 displays how the evolution of entanglement 
depends on the dipole-dipole strength for n = 1. From these figures, we find that 
entanglement revealed by the solid lines (a = 0.1) is stronger than that by the 
dotted lines (a = 1) at most time. The value of atomatom negativity is larger 
when the dipole interaction between the atoms increases. This dependence is 
much sharp for incoherent initial atomic state. 

4. Conclusion 

In this paper we have investigated the effect of the atomic coherence on 
the entanglement of two dipole coupled two-level atoms when only one atom is 
trapped in a lossless cavity and interacts with one-mode thermal field, and the 
other one can be spatially moved freely outside the cavity. We have shown that 
thermal field can produce atom-atom entanglement for all pure initial atomic 
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states. The results also show that the atom-atom entanglement can be con¬ 
trolled by changing the system parameters, such as the amplitudes of the polar¬ 
ized atoms, the mean photon numbers of thermal field, and the strength of the 
dipole interaction (or distance between atoms). We have derived that for con¬ 
sidered model the dipole-dipole interaction can produce the appreciable amount 
of entanglement only for small thermal noise if the atoms are prepared in the 
coherent states. 
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Appendix: Expressions for density matrix Elements in Eq. (11) 

The evident expressions for matrix elements in (11) are 


/Oil = H 2 X! P "I^41,™| 2 + H 2 ^2Pn\Z42,n-l\ 2 + be* ^ Pn-^42,n-l-^43,n_l + 

n —0 n=1 n=l 

+cb* EE Pn^43,n— 1 ^42,n— 1 1^1 ^ ^ Pn \ ^43, ti— 1 1 H - |^| ^ ^ Pn | ^44, ti—2 1 •> 

71=1 71=1 71=2 

Pl2 — ^ ^ Pn^41, 71^22,71—1 ^ ^ Pn ^41,n ^23,n— l”l” 

71=1 71=1 

+bd* '^2 PnZ42 ,n—l Z24,n—2 + C(f p n Z43 t n-lZ24 n _2 + 
n —2 n =2 

+Pi{bd*Z4 2 i oG 24 + cd*Z 43 qG 24 ) + Po{ab*Z 41 qG 22 + ac*Z 41 ^ 23 ), 

Pl3 = (lb ^ ' Pn-^41,n-^32,n— 1 ~i~ CtC ^ ' Pn-^41,n^33,n—l~f~ 

71=1 71=1 

+bd* p n Z42 t n-lZ^_ n _2 + cd* Pn-^43,ra-1^34. n - 2 + 

ra=2 n —2 

+Pi(bd* Z 42 fiG ^ 4 + cd* Z43 ,qG 3^) + po(ab* Z410G32 + dc* Z4ipG^ 3 ), 

p \4 = ad’ *E 2^41,71^14,71-2 H- Pi o-d* ^4i,iG^4 + pqclcI* 

n —2 
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P22 = |«| 2 y^Pn|-^21,n| 2 + M| 2 ^ Pn|-^24,ra-2|~ + 
n=0 n—2 

+0C* Pn-^22,n—1-^23,n—1 + c ^* Pn-^23,71-1-^22,71-1 + 

n=l n=l 


+ |c| 2 Pn|-^23,71-1 | 2 + |0| 2 Pn\^22,n-1 1 2 + 

n—1 7i=l 

+Pi|rf| 2 |G 2 4| 2 +Po(|6| 2 |G 2 2| 2 + 6c*G22G*3)+po(c6' t G23G* 2 |c| 2 |G23| 2 ), 

023 = |«| 2 y; PnZ2\,nZ^i n + \d\ 2 p n Z24,n-2Z^ n _ 2 + 

n=0 n— 2 

+00* y p„Z22,w-l-^33,n-l + 00* ^ Pn^23,ra-1^32, n -l + 

n=l 7i=l 

+w 2 E Pn-^23,n-l-^33,n-l + |0|~ y Prt^22,ra-1^32.n-l + 

n=1 n= 1 

+p[ 1]M| 2 G 2 4G^3 + p[O](|0| 2 G 22 G^ + bc*G 22 G* 3 3 + cb*G 23 G* 32 + |c| 2 G 23 Gy, 

033 = H 2 y Pn|^31,nl 2 + |0| 2 E Pn |^32,n-l 1“ + 00* y P n ^32,ra-1^33,n-l + 
n—0 n=l ti =1 

+c0* y PnZ 33,71—1-^32,71—1 l C l 2 E Pn|-^33,n-l | 2 2 + |d| 2 y Pn|-^34,n-2| 2 + 

n=l n—1 ti=2 

+PiM| 2 |G34| 2 + po(|0| 2 |G 32 |“ + be* G32G33 + c 6*G 3 3G 2 2 + | c| 2 1G331 2 , 

034 = ab* 'EPnZsi.nZu.n-l + CLC* Pw^31,n^l3,rt-1 + 
n—1 n=l 

+ 00* yi Pn-^32,n-l-^14,n-2 + 00* El Pn-^33,n-l-^14,n-2 + 

71=2 71=2 

+Pi(M*Z 32> oGl4 + cd*Z3 3> oG* 4 ) + po(a0*Z 3 i jO G; 2 + 00*^31^0*3)+ 
+Po(00*G 32 + cd*G 33 , 

044 = \a\ 2 y Pnl-^llj 2 + |^| 2 ypnl^-ll 2 + be* PnZ\ 2 , n-l-^13,rt-l + 

71=0 71=1 71=1 

+00* y p n Zi 3) n-lZl 2n _i + |c| 2 y Pn\Zl3,n-l\ 2 + |0| 2 Pn | ■Z’l4,n-2 | 2 + 

71=1 71=1 71=2 

+P 1 MIIG 14 I 2 + Po(|^| 2 |Gi 2 |“ + 6c*Gi 2 G* 3 + c6*Gi 3 G* 2 + |c| 2 IG 13 1 2 + |0| 2 ), 


where 

0 = 008 01008 02, b = cos 0i sin0 2 e* V2 , c = cos0 2 sin0ie* Vl , 0 = sin0i sin0 2 e^ Vl+V2 ^. 
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Figure 2: The negativity as a function of gt for the model with n = 1 (solid) and n = 2 
(dashed). The initial atomic state: a) |^(0 ))a = |+, —) (a) and |^(0))i = (l/\/2(|+) + 
|—)), |^(0 ))2 = (l/-\/2(|+) + |—)) (b). The dipole strength a = 0.5 
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Figure 3: The negativity as a function of gt for the model with a = 0.1 (solid) and a = 1 
(dashed). The initial atomic state: a) |\I/(0))^ = |+, —) (a) and |T(0))i = (l/\/2(|+) + 
|—)), |^(0 ))2 = (l/\/2(|+) + |—)) (b). The mean photon number n = 0.5 
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